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Status Report 
1 August 1977 

This status report on "Low Cost Avionics for the Global Positioning System," 
NASA Grant NSG 1374, is forwarded at this time to report conclusions on several 
recent investigations Also, ongoing work is detailed. 

Costas Loop and Correlation Analyses 

The key to feasibility and usefulness of the single-channel GPS receiver is 
the speed with which code and carrier synchronization may be established as the 
receiver cycles over a sequence of satellite signals. Thus a considerable 
amount of effort has been spent on modeling and simulating the Costas Loop, 
together with a fast algorithm for estimating and forcing "lock". 

Enclosed is a copy of TCSL Research Memorandum 77-06, "Costas Loop and 
Correlation Analyses for the Low Cost GPS Receiver." This report is composed 
of two engineering reports prepared by Captain David A. McClung, U.S. Army, 
in partial fulfillment of the requirements for the degree. Master of Engineering, 
at Texas A&M University The report covers the mathematical modelling and 
Monte Carlo (FORTRAN) simulation of two types of Costas Loops in noise, including 
a correlation detector and algorithms for sweeping the loops into lock. 

The initial approach to the acquisition problem, which McClung has followed, 
is to sweep the Costas Loop VCO as rapidly as possible over a 13 Khz. fre- 
quency uncertainty range, and estimate the frequency at which maximum correlation 
occurred Then the VCO is reset to the estimated frequency and the loop is 
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<ill()W('d to "()ull" into lock naturally. In the limited time available to him 
he was able to obtain simulation results only for one particular sweep rate 
and one particular value of correlator bandwidth. 

McClung's first results show that a 10 Khz. uncertainty range can be 
swept in 70 milliseconds However, his noiseless linear estimate of maximum 
correlation frequency had a worst case error of 150 Hz. With this tuning 
error, the loop pull-in time is 12 and .42 seconds for a 100 Hz. loop band- 
width, for the Sine Loop and Tangent Loop respectively With a 70 milli- 
second sweep and 2048 trial code correlations, the maximum time aevoted to 
sweeping would be 141 seconds. So 0.5 seconds for the loop to pull in to 
final lock does not seem comparatively excessive. However, the total worst 
case time to lock of 141.5 seconds is absolutely excessive for nominal use. 

The design goal is for cycle times between satellites on the order of 1 second. 

It is apparent that in order to reach the design goal of 1 second 
between satellites for nominal operations, both the frequency uncertainty and 
the code delay uncertainty must be drastically reduced. Only by limiting the 
range of the frequency sweep and the number of trial correlations can the desired 
acquisition time be realized To reduce the sweep range and number of trial 
correlations will require computational aid from the receiver's microcomputer 
processor 

After each satellite has been acquired once, the 50 BPS data read, and the 
satellite ephemendes stored, then the frequency uncertainty due to satellite 
doppler may be largely eliminated. This requires that a rough estimate of the 
users position be available At that point the A Prion frequency uncertainty, 
or sweep range, of 13,532 Hz can be reduced by 9767 Hz. to 3,. ’65 Hz. Since 
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a relative frequency estimator is used in the Costas Loop it is quite likely 
that the sweep range can be further reduced by partially eliminating the un- 
certainty due to oscillator long term stability (945 Hz.) and aircraft doppler 
(2820 Hz.) 

Assuming a one-second cycle time between satellites, a sequence of four 
satellites and a worst case aircraft velocity of 600 statute miles per hour, 
the worst case code delay uncertainty for two consecutive acquisitions of the 
same satellite is 2144 meters or 7.154 microseconds, without knowledge of the 
aircraft velocity vector. This is equivalent to 7.32 PN bits or 14.6 trial 
correlations 

Suppose that, conservatively, we set the code delay uncertainty at 20 
trial correlations and the frequency uncertainty at 3,765 Hz. (the satellite 
Doppler uncertainty has been eliminated). Let us allow 0.5 seconds for 20 
trial correlations and frequency sweeps. This requires a sweep time of 25 
milliseconds or 150,000 Hz/sec sweep rate, which is approximately the rate 
demonstrated by McClung Thus, McClung's results are applicable to the single- 
channel GPS receiver under realistic assumptions for Doppler/Deluy-aiding 
acquisition 

The work which will be pursued in the Costas Loop area during the remainder 
of Phase I (until October 1, 1977) will be to refine the lock- frequency esti- 
mation algorithm and to evaluate its performance in the noise environment to 
be expected under nominal conditions, which is for a signal to noise spectral 
density of 38 4 dB Hz The correlation filter bandwidth may be reduced somewhat 
from the value demonstrated by McClung, 314 Hz , depending on performance 
Also, the Costas-Loop bandwidth may be adapted from its acquisition value of 
100 Hz to Its tracking value of 10 Hz. during the maximum pjll-in time. It 
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will be determined whether to acquire the loop in the Sine configuration or in 
the Tangent configuration In addition, the hold-in behavior of both loops 
will be evaluated when they are subjected to phase accelerations (frequency 
ramps) after being initially locked. This evaluation will determine whether the 
Tangent Loop is more immune to dynamics than the Sine Loop. 

Keplenan Orbit Simulation and the Navigation Algorithm 

Ephemendes data for the 24 GPS satellites had previously b'^en obtained 
from a circular deterministic orbit algorithm obtained from Dr. Thomas Wu and 
Major Kenneth A. Myers of the USAF Avionics Laboratory at WPAFB, Ohio. Actual 
satellite data transmitted in Data Block II are the Keplenan orbit parameters 
that are required to solve the Keplenan orbit equations in an ECEF (earth 
centered earth fixed) coordinate system 

The satellite orbit parameters range from 8 to 32 bits of accuracy. 

Memory constraints and computation time dictated a desire to create a micro- 
processor architecture based on a 32 bit floating point system (24 bit mantissa 
and 8 bit exponent and sign) Thus a major purpose of this research was to 
determine the affect of truncating 32 bit data parameters down to 24 bits. 

The initial simulation produced ephemens data for a single satellite 
in one complete 12 hour orbit Two runs were made one with 32-bit precision 
and one with 24-bit precision. The RMS difference was defined as the distance 
error of the 24-bit implementation. The maximum error encountered in a 12-hour 
orbit IS 4990 feet. 

A subsequent simulation was performed to determine the smoothing effect 
(if any) of the Navigation Algorithm on the ephemens error. A simulated USAF 
C5A aircraft flight from California to Hawaii showed fixes taken at 30 second 
intervals to be less than 100 feet in error 82% of the time and less than 325 
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feet all of the time. These results compare almost exactly witf the previously 
produced noisy model simulation for the circular orbit algorithms 

These findings suggest that the noise completely obscures the truncation 
error of the algorithms, however, we have subsequently determined that the 
reverse is actually the case. 

Double precision (56 bit mantissa) runs were made without noise substan- 
tiating algorithm accuracy obtained at WPAFB previously. Single precision 
(24 bit mantissa) runs were made without noise and nearly identical results 
were obtained as the runs made with noise That is, typical errirs of 100 feet 
or less occur in Q2% of the fixes with or without noise if only 24 bits of 
precision are used in the computations. Hence, it is conclusive that (with 
the noise models used) truncation (or round-off error) is the dominant factor 
of error production for a 24-bit computer system in the low cost GPS receiver. 

The existing navigation software has just been converted to operate on 
the DECll/40 USAF Computer at Texas A&M. Experiments have been defined that 
will be used to modify the navigation software from its present parallel 
input of satellite information to a low-cost single-channel serial input 
format. Satellite selection will also be studied The existing algorithm for 
selection will be used to obtain a set of four satellites with good GDOP, 

Then these four satellites will be retained until one (or more) disappears 
from view GDOP and position error will be recorded to determine the effective- 
ness of this procedure A new selection procedure will be used when a satellite 
disappears to determine if any of the original set can be re ained. Again 
GDOP and position error will determine the effectiveness of the procedure. 



Hardware Effort 


The hardware research program has concentrated on the microprocessor- 
controlled phase shift circuitry that permits the locally generated Gold 
code to be cross-correlated with the received Gold code in order to obtain 
range. The phase shift circuitry includes a divide-by-44 counter, a digital 
pulse cancellation subsystem, and a variable-phase divide-ty-35 counter 
system These three systems have been constructed in prototype form by 
Mr 0 Ladipo, an Electrical Engineering graduate student. Typical 7400 
TTL digital devices were used, limiting the prototype to several mHz in 
frequency The performance was satisfactory at that frequency range, 
however, and translation to the 35 mHz range is not expected to be difficult. 
Mr. Ladipo has completed an engineering report on his work, and copies 
of that report are enclosed 

Additional study has been directed toward the generation of the various 
Gold codes used by each of the 24 satellites. In the next reporting 
period the generating polynomials for each code will be tested, and 
hardware/software subsystems for generating these codes under micro- 
processor control will be developed Also, a systematic approach to 
initializing the sequential code generator to any desired phase will 
be developed 
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Selection of Receiver Parameters 


Basic Assumptions 

Worst Case Satellite Radial Velocity 2,078 MPH (stat.) 

Highest Aircraft Velocity 600 MPH (stat.) 

Oscillator Long Term Stability (Aircraft): 3 x 10"^ 


Calculations 

The oscillator output frequency is 

= F + A - € c=3x 10’^ X F (1) 

where 

F = = 35805000. Hz. (2) 

and A IS chosen as a design parameter. Note that 44F is the satellite trans- 
mitter frequency, equal to 1575.42 Mhz. 

The second I.F frequency, which is the Costas Loop input frequency is 
chosen to be less than or equal to 30 Khz. under all conditions. This is to 
give the possibility of eventual total digital realization for the Costas Loop 

The normalized radial velocity between aircraft and satellite (worst case) is 


V 

c 



Vs 600 + 2078 
c ■ 186219 X 3600 


= 3 99 X 10'® 


The total tuning increment for the Costas Loop is given by 

Af = (2)(44)(e + F(^)) 

= 88(3 X 10"^ + 3.99 x 10"®)35805000 
= 13532 Hz 


(3) 


(4) 
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The lower end of the second I.F band is then 


30,000 - 13,532 = 16,468 Hz 
and the center of the band is 

(30,000 + 16,468) - 2 = 23,234 Hz 
Thus A IS solved for as 

A = = 528 0 Hz. (5) 

The first I F passband is defined by 

f^^ = F + 45a - (45e + 44 F(^)) 

= 35828760 - 6,769 (6) 

Thus, the I F center frequency and bandwidth are respectively 

F^^ = 35.828760 Mhz 
BW = 13 5386 Khz 

The relationship between first I.F bandpass and oscillator f req lency is shown 
in Figure 1 


1 4 17,001 4 13,538 — ^ 



35.80500 35.828760 


Figure 1. Relation of I F. to Oscillator Frequencies 
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F+45A - (45e + 44F(^)) 

35.828770 Mhz 
Bw = 13.6 Khz. 


44(A t (e + F(^)) 
16,460 - 30,000 Hz. 



F = 35.805000 Mhz. 
v/c = 3.99 X 10"® 
F(v/c) = 1^3 Hz. 

A = 528.0 Hz. 
e = 10.7 Hz. 


I G PS RecelvefI 
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LINK MARGIN ANALYSIS 

1. Transmitter Power 20 Watts (Radiated) 

2. Transmitter Circuit Loss 

3. Transmitter Antenna Gain 

4. Transmitter Antenna Pointing Loss 

5. Space Loss f - 1575.4 Mhz., r = 25824 km. 

6. Antenna Polarization Loss 2 db ellipticity XMT 

7. Receiving Antenna Gam 3 db ellipticity Rev. 

8. Receiving Artenna Pointing Loss 

9. Receiving Ci.'cuit Loss 

10. Net Loss Link (Add 2 through 9) 

11. Total Received Power (Add 1 plus 10) 

12. Receiver Noise Spectral Density T = 625 

F = 5 db 

13. Received Signul/NSD (11 minus 12) 

COSTAS LOOP 

14. Costas Loop BW 10 Hz. 

15. Costas Loop SNR (13 minus 14) 

16. Required Loop SNR 

17. Carrier Margir (15 minus 16) 

DATA 

18. Bit Period sec. 

19. Received Energy/NSD, E/No (13 plus 18) 

20. Required E/No 

21. Data Margin C9 minus 20) 
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